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ABSTRACT 

A  comparison  of  two  different  autopilot  configurations  and  their  effect  on  missile 
response  is  presented.  The  comparison  includes  the  steps  taken  in  determining  missile 
parameters  from  wind  tunnel  data  and  flight  condition  data.  The  missile  parameters  are 
coupled  with  two  different  autopilot  configurations  to  determine  any  significant  advan- 
tage of  one  configuration  over  the  other.  Pole  placement  is  used  in  determining  required 
autopilot  feedback  and  feed  forward  gains.  Simulations  of  each  autopilot  are  conducted 
and  the  responses  are  compared. 
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I.     INTRODUCTION 


The  purpose  of  an  autopilot  in  a  missile  is  to  cause  the  missile  to  maintain  stability 
as  it  travels  along  its  flight  path.  Commonly  the  autopilot  provides  stability  about  the 
missile's  roll,  yaw  and  pitch  axes.  Autopilots  generally  use  displacement  gyros  with 
feedback  being  provided  by  electronic  sensors  which  detect  axial  displacement  due  to 
missile  rotation  about  one  or  more  of  the  above  mentioned  axes. 

Autopilots,  as  with  most  other  man-made  devices,  can  be  constructed  in  a  variety 
of  diflcrent  configurations.  Some  configurations  may  have  advantages  over  others  in 
terms  of  parameters  such  as  cost,  weight  or  effectiveness,  to  name  a  few. 

In  this  research  paper  two  pitch  axis  autopilot  configurations  are  compared.  The 
first  is  a  three  loop  autopilot  configuration,  employing  body  rate  feedback,  lateral  ac- 
celeration feedback  and  an  additional  synthetic  stability  feedback  loop  as  depicted  in 
Fie ure   1. 
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Figure   1.      Autopilot  Configuration  1 


The  second  is  a  two-loop  autopilot  configuration,  employing  proportional  plus  integral 
compensation  in  the  feed  forward  path  with  bod}'  rate  and  lateral  acceleration  feedback 
as  shown  in  Fieure  2. 
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Figure  2.      Autopilot  Configuration  2 


In  these  two  figures  >/c  is  the  commanded  lateral  acceleration  in  the  pitch  plane,  ;/iM 
is  the  measured  lateral  acceleration  in  the  pitch  plane,  S  is  the  missile  control  surface 
dcllection  angle  and  0  is  the  missile  pitch  angle  rate.  All  the  subscripted  K  variables  are 
gains  which  arc  constant  at  each  design  point. 

The  purpose  of  the  comparison  of  the  two  autopilot  configurations  is  to  determine 
what  advantage,  if  any,  one  may  have  over  the  other  in  terms  of  missile  response. 
Comparisons  arc  made  of  the  two  autopilots  for  three  different  sets  of  flight  conditions. 
Cach  set  of  flight  conditions  is  referred  to  as  a  design  point.  For  example,  design  point 
one  will  refer  to  the  first  set  of  flight  conditions.    These  flight  conditions,  or  design 


points,  are  specified  in  Chapter  Two.  The  desired  autopilot  response  for  each  config- 
uration is  given  as  approximately  10  radians  second  with  a  damping  coefficient  of  0.5. 

Pertinent  missile  parameters  are  derived  in  Chapter  Two.  using  given  wind  tunnel 
data  and  flight  condition  data.  Any  needed  parameters  which  cannot  be  derived  or 
otherwise  obtained  through  availible  information,  are  assumed  using  sound  engineering 
judgment. 

In  Chapter  Three  the  missile  and  autopilot  transfer  functions  are  developed.  The 
feedback  gains  necessary  to  meet  the  desired  autopilot  design  specifications  are  obtained 
allowing  development  of  simulation  studies. 

Simulations  are  developed  and  conducted  using  Dynamic  Simulation  Language 
(DSL)  in  Chapter  Four  and  the  results  are  analyzed. 


II.     DEVELOPMENT  OF  MISSILE  PARAMETERS 


A.     INTRODUCTION 

In  order  to  compare  the  effects  that  the  two  autopilot  configurations  have  on  missile 
performance,  determination  of  certain  missile  parameters  must  first  be  made.  Specif- 
ically, the  transfer  functions  0(s)jd(s)  and  i]M(s)jS(s),  which  are  the  missile  dynamics 
blocks  in  the  autopilot  block  diagrams  of  Figure  1  on  page  1  and  Figure  2  on  page  2, 
must  be  obtained. 

In  order  to  determine  these  transfer  functions,  wind  tunnel  data  and  three  sets  of 
flight  condition  data  arc  provided  later  in  this  chapter.  The  procedures  for  processing 
this  data  are  also  presented  later  in  this  chapter. 

Prior  to  developing  the  missile  parameters,  an  explanation  of  the  notation,  termi- 
nology, and  coordinate  system  used  is  in  order.   [Ref.  I:  p.  132]  Refer  to  Figure  3. 


Figure  3. 


Missile  Notation 


The  notation  in  Figure  3  represents  the  following: 

CG  is  the  center  of  gravity  of  the  missile. 

CP  is  the  aerodynamic  center  of  pressure  of  the  missile. 

VM  is  the  velocity  of  the  missile. 

SRL  is  the  spatial  reference  line,  a  reference  fixed  in  space. 

FL  is  the  lift  force.    It  is  perpendicular  to  the  velocity  vector. 

W  is  the  weight  of  the  missile. 

MP  is  the  pitching  moment.    It  is  taken  as  positive  in  the  counter-clockwise  direction. 

a  is  the  angle  of  attack.    It  is  the  angle  taken  from  the  velocity  vector  to  the  longi- 
tudinal missile  axis. 

6  is  the  angle  between  the  SRL  and  the  longitudinal  missile  axis. 

S  is  the  fin  deflection  angle.  It  is  taken  from  the  longitudinal  missile  axis  to  the  control 
fin  axis. 

ft  is  the  angle  between  the  gravity  vector  and  the  SRL. 

Other  notation  and  terminology  will  be  discussed  as  it  is  introduced.  It  is  important  to 
remember  that  this  problem  considers  only  motion  in  the  pitch  plane  but  that  similar 
methods  may  also  apply  to  yaw  motion. 

B.     FLIGHT  CONDITION  DATA 

The  three  sets  of  flight  condition  data  given  in  Table  1  on  page  6  establish  the  de- 
sign points  around  which  the  autopiloted  missile  will  be  examined.  The  flight  condition 
data  will  be  used,  in  conjunction  with  wind  tunnel  data  presented  in  the  next  section,  to 
establish  the  missile  parameters  and  autopilot  gains  necessary  to  develop  and  conduct 
simulations  of  the  autopiloted  missile. 

Some  manipulation  of  the  flight  condition  data  is  required  in  order  for  it  to  be  used. 
These  conversions  and  calculations  will  be  made  as  necessary  and  explained  when  used. 


Table  1.     FLIGHT  CONDITIONS  AT  THREE  DE- 
SIGN POINTS 


DESIGN  POINT  I 

T 

= 

Time  (referenced  to  launch)  =  9.5  s 

Alt 

= 

Altitude  : 

=  4.993  km 

V 

= 

Velocity  = 

=   1112.4  ms 

M 

=  Mass  = 

450  kg 

q 

= 

Dynamic 

Pressure  =  455642  N/m2 

I 

= 

Moment  o 

f  Inertia  =  720  kg  m2 

Powei 

On 

CP 

is 

265  cm  aft 

of  nose 

DESIGN  POINT  2 

T 

= 

18.5  s 

Alt 

= 

13.54  km 

V 

= 

1086.4  m 

s 

M 

=  370  ks 

q 

= 

146829  N 

m2 

i 

= 

693  k2  m2 

Powei 

■Off  * 

CP 

is 

265  cm  aft 

of  nose 

DESIGN  POINT  3 

T 

= 

22.5  s 

Alt 

= 

17.S2  km 

V 

= 

934.5  m.'s 

M 

■■  365  ks 

q 

= 

54637  N/rrf 

i 

= 

6S7  kg  m: 

Powei 

■  Oil 

CP 

is 

266  cm  aft 

of  nose 

C.     WIND  TUNNEL  DATA 

The  wind  tunnel  data  presented  in  this  section  will  be  instrumental  in  determining 
missile  parameters.  The  data,  shown  in  Table  10  on  page  37  through  Table  13  on  page 
40,  are  extracted  from  the  complete  wind  tunnel  data  set  for  the  missile  used  in  this  re- 
search paper.  These  data,  as  well  as  the  three  sets  of  flight  condition  data  contained  in 
Table  1,  were  provided  by  the  U.S.  Army's  Missile  and  Space  Intelligence  Center 
(MSIC). 

The  first  two  sets  of  wind  tunnel  data  are  the  axial  force  coefficients.    These  are 
given  in  Table  10  and  Table  11  for  the  Power  On  and  Power  Off  cases. 


The  third  set  of  wind  tunnel  data  is  the  normal  force  coefficients.  These  data  in 
Table  12  on  page  39  will  be  used,  along  with  the  axial  force  coefficient  data  to  calculate 
the  lift  and  drag  of  the  missile  at  the  design  points. 

The  final  set  of  wind  tunnel  data  is  the  pitching  moment  coefficients.  These  data. 
in  Table  13  on  page  40  will  be  used  to  determine  the  pitching  dynamics  of  the  missile. 

D.     WIND  TUNNEL  DATA  INTERPOLATION 

To  enter  the  wind  tunnel  data  at  the  desired  design  point,  the  Mach  number  must 
be  calculated  from  the  flight  condition  data.  Using  given  missile  velocity  and  altitude, 
with  temperature  obtained  from  gas  tables  and  graphs  found  in  [Refs.  2,3],  the  Mach 
number  for  each  design  point  is  calculated  by  gas  law  equation  techniques.  These  cal- 
culations are  not  included  here.  The  Mach  number  at  each  design  point  is  given  in 
Table  2. 

Table  2.     MACH  NUMBERS  AT  DESIGN  POINTS 


DESIGN  POINT 

MACH  NUMBER 

1 

3.48 

j 

3.68 

■^ 

3.17 

The  availability  of  the  Mach  number  at  each  design  point  allows  entry  of  the  wind 
tunnel  data  and  interpolation  of  the  various  coefficients  for  each  design  point.  The  in- 
terpolation of  the  wind  tunnel  data  at  each  design  point  is  accomplished  using  a  simple 
linear  interpolation  scheme.  The  interpolation  provides  the  coefficients  at  each  design 
point.  These  interpolated  coefficients  are  given  in  Table  14  on  page  41  through 
Table  17  on  page  43. 

E.     CALCULATION  OF  LIFT  AND  DRAG  COEFFICIENTS 

Once  the  axial  and  normal  force  coefficients  at  the  design  points  have  been  deter- 
mined the  lift  and  drag  coefficients  can  be  calculated.  The  lift  force,  FL,  and  drag  force, 
FD  ,  are  related  to  the  normal  force,  FN,  and  the  axial  force,    FA,    by  the  trigonometric 


relationships 


FL  —  F^-cosa.  —  FAs'ma 


(2.1) 


and 


FD  =  Fss'ma  +  FAcosa. 


(2.2) 


These  relationships  are  depicted  in  Figure  4  below. 
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Figure  4.      Force  Relationships 

Since  the  cocITicients  are  dimensionless  representations  of  the  forces,  the  same  re- 
lationships apply.    Specifically, 


C^  =  CA<:osa  —  C^sina 


(2.3) 


and 


C^  =  C/Vsina  +  C^cosoc 


(2.4) 


where  CL  is  the  lift  coefficient,  CD  is  the  drag  coefficient,  CA  is  the  axial  force  coefficient, 
and  Cv  is  the  normal  force  coefficient.   [Ref.  1:  p.  134] 

The  application  of  these  equations  to  the  interpolated  wind  tunnel  data  in  Appendix 
B  renders  lift  and  drag  coefficients  at  the  design  points.  These  lift  and  drag  coefficients 
are  tabulated  in  Appendix  C. 

F.     ANALYSIS  OF  WIND  TUNNEL  DATA 

There  are  four  missile  parameters  which  are  of  interest  in  the  development  of  the 
missile  dynamics  transfer  functions.  These  parameters  are:  Fu,  which  is  the  lift  force  per 
angle  of  attack,  Fu,  which  is  the  lift  force  per  fin  deflection  angle,  A/Pa,  which  is  the 
moment  about  the  center  of  gravity  per  angle  of  attack,  and  MpS,  which  is  the  moment 
about  the  center  of  gravity  per  fin  deflection  angle. 

These  parameters  can  be  determined  from  their  corresponding  coefficients  and  a  few 
physical  dimensions  of  the  missile.  The  corresponding  coefficients  are  discussed  in  the 
paragraphs  which  follow. 

CLx  is  the  lift  coefficient  per  angle  of  attack.  It  is  determined  by  graphing  the  lift 
coefficients  against  the  angle  of  attack  and  approximating  the  slope  of  the  nearly  linear 
graph. 

Clf  is  the  lift  coefficient  per  fin  deflection  angle.  It  is  determined  by  graphing  the  lift 
coefficients  against  the  fin  deflection  angle  and  approximating  the  slope. 

Similarly,  CM,  and  Cm  are  the  pitching  moment  coefficients  per  angle  of  attack  and 
fin  deflection  angle,  respectively.  They  too  are  obtained  by  graphing  the  coefficients 
against  the  angles  and  approximating  the  slope  of  the  resulting  graph. 

The  graphs  of  all  the  coefficients  versus  corresponding  angles  are  contained  in  Ap- 
pendix D.  The  graphs  are  made  for  each  set  of  coefficients  at  each  of  the  three  design 
points.    The  coefficients  CLa,  CLi,  CM%  ,  and  Cm  are  tabulated  in  Table  3  on  page  10. 


Table  3.     COEFFICIENTS  DUE  TO  ANGLE  OF  ATTACK  AND 
FIN  DEFLECTION  ANCLE 


DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 
POINT  3 

Ci* 

0.27 

0.25 

0.28 

Cu 

0.03 

0.02 

0.03 

Cms 

-0.15 

-0.12 

-0.20 

Cm 

-0.17 

-0.15 

-0.22 

G.     CALCULATION  OF  MISSILE  PARAMETERS 

As  mentioned  in  the  preceding  section,  the  parameters  of  interest  are  Fu,  FLI,  MPr, 
and  MFi.  With  knowledge  of  the  lift  and  moment  coefficients  previously  obtained,  and 
knowledge  of  some  basic  physical  dimensions  of  the  missile,  these  parameters  are  deter- 
mined. The  necessary  missile  dimensions  are  shown  in  Figure  5.  These  missile  dimen- 
sions were  provided  by  the  Missile  and  Space  Intelligence  Center. 


SIDE  VIEW 


5.3  m 


.DIAMETER  0.38  m 


END  VIEW 


Figure  5. 


Missile  Dimensions 
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The  equations  used  to  determine  the  missile  parameters  from  the  coefficients  are: 

FLoc=CLo:CIA'  i2-5) 

FL6=CL6qA,  (2.6) 

mp*=cm*1Ac>  (2-7) 

and 

MPS=CmqAc  (2.8) 

where  q  is  the  dynamic  pressure,  provided  in  the  flight  conditions.  A  is  the  missile  char- 
acteristic area,  and  c  is  the  characteristic  length.   The  characteristic  area,  A,  is  taken  as 

A  =  n^-  (2.9) 

where  d  is  the  missile  diameter.  38  cm,  as  given  in  Figure  19  above.  The  characteristic 
length  is  taken  as  the  missile  diameter.   [Ref.  1:  p.  134] 

Applying  the  above  equations  to  the  known  coefficients  produces  the  missile  pa- 
rameters at  the  three  design  points.  Use  of  these  parameters  in  later  chapters  requires 
that  they  be  given  in  force  or  moment  per  radian  rather  than  per  degree.  This  conver- 
sion is  made  mulitplying  each  parameter  by  ISO/^  .  These  converted  parameters  are 
tabulated  in  Table  4  on  page  12. 

Knowledge  of  these  missile  parameters  will  enable  the  derivation  of  the  missile  dy- 
namics transfer  functions  in  the  following  chapter. 
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Table  4.     MISSILE  PARAMETERS  AT  DESIGN  POINTS 


PARAMETER 

DESIGN 
POINT  1 

DESIGN 

POINT  2 

DESIGN 

POINT  3 

Fu 

799276  N 

238350  N 

99695  N 

FL6 

88808  N 

19080  N 

10657  N 

Afft 

-168736  Nm 

-43545  Nm 

-26986  Nm 

MP6 

-191368  Nm 

-54431  Nm 

-29679  Nm 
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III.     DETERMINATION  OF  TRANSFER  FUNCTIONS 

A.  OVERVIEW 

In  this  chapter,  the  missile  parameters  obtained  in  the  previous  chapter  are  applied 
to  the  equations  of  motion  which  govern  the  flight  of  the  missile.  From  this,  the  missile 
dynamics  transfer  functions.  6(s)IS(s)  and  riM(s)ld(s),  are  derived.  By  then  inserting  these 
two  missile  dynamics  transfer  functions  into  the  two  autopilot  transfer  functions,  the 
complete  autopilot  transfer  functions  can  be  expressed  in  terms  of  missile  parameters 
and  feedback  gains. 

B.  MISSILE  DYNAMICS  TRANSFER  FUNCTIONS 

Consider  the  forces  and  moments  acting  on  a  missile  while  in  flight.  These  are  de- 
picted in  Figure  3  on  page  4.  Summing  the  forces  in  the  transverse  direction,  that  is 
perpendicular  to  the  velocity  vector,  yields 

Yft=  ih!m  =  FL-  IV sin(/?  +  6  -  a).  (3.1) 

Since 

FL  =  FLS.  +  FL6(S  +  z),  (3.2) 

the  substitution  is  made  yielding 

*lMm  =  FuP-  +  F^S  +  a)  -  Wsin(P  +  6  -  a).  (3.3) 

Dividing  through  by  the  mass,  m,  yields  the  lateral  acceleration,  y\M  . 


>?A,  =  -^a+-^(<5  +  «)--7FSin(£  +  0-a)  t3'4) 


The  weight  of  the  missile  is  neglected  at  this  point  for  simplicity.  The  missile's 
weight  would  affect  the  lateral  acceleration  of  the  missile  most  when  the  angle 
(/?  +  6  —  a)  is  90  or  270  degrees.  At  this  angle  the  contribution  to  lateral  acceleration 
would  be  plus  or  minus  1  G.    This  effect  would  be  offset  by  fin  deflection  and  angle  of 
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attack  trim  angles  of  much  less  than  one  degree  at  any  of  the  three  design  points. 
Neglecting  missile  weight  leaves  the  approximation 

im-  -ura+-zr(d  +  (/-)  =  c'- m +  d  ~ur-  ( 3-^ 

which  is  used  to  determine  rjM  at  each  design  point  in  terms  of  a  and  S. 

In  a  similar  manner,  summing  moments  about  the  center  of  gravity  of  the  missile 
with  the  counter-clockwise  direction  being  referenced  as  the  positive  direction  yields 

Ymcg  =  16  =  MP  (3.6) 

where  /  is  the  missile's  moment  of  inertia,  6  is  the  missile's  angular  acceleration  in  the 
pitch  plane,  and  MP  is  the  pitching  moment.  The  pitching  moment,  MP,  is  given  by  the 
equation 

MP  =  aAfrt  +  6Mpe  +  aAfft,  +  (a  +  S)MP6.  (3.7) 

Since  a.  ^  6  and  a  =s  6  the  substitution  is  made  leaving 

MP  =  «(MPi  +  MJ  +  a(MPa  +  MP8)  +  S(MP6).  (3.8) 

Combining  and  rearranging  the  above  equations  yields 

..      W>  +  Up*)   .       (A/fa  +  Mm)  MP6 

a— a— - a  = — - — o.  (3.9) 

In  Equation  (3.9),  the  term 

(MPi  +  Mpb) 


is  the  aerodynamic  damping  term.  It  is  beyond  the  scope  of  this  research  paper  to 
properly  obtain  this  term.  It  is,  therefore,  assumed  that  the  missile  is  designed  with  a 
damping  coefficient  of  £  =  0.5.  The  aerodynamic  damping  term  will  be  different  nu- 
merically at  each  of  the  three  design  points  and  will  be  referred  to  hereafter  as  D.   It  will 
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be  calculated  later  in  this  chapter.   [Ref.  4] 
Substituting  D  into  Equation  (3.9)  yields 

a  —  Da.  — a  = — - — <5.  (3.10) 

Using  LaPlace  transforms  and  assuming  zero  initial  conditions  this  equation  transforms 
into 

5  a(s)  -  Z)sa(s) u(s)  =  — - — 5(s).  (3.11) 

Rearranging,  this  becomes 

M 


P6 

cr.(s)  r 


S{s)         2                (M^  +  Mw) 
5    —  L>5  — — : 


(3.12) 


Multiplying  both  sides  by  5  renders 


MP6 


sois)  I  vis) 


d{s)  2  {MPx  +  MM)        5{s) 

S    —  Ds  — : 


(3.13) 


Since  a.  =  6, 


vis) 

6(s) 
'  5{s)  - 

s 

MP6 
I 

d(s)  ' 

0 

5"  ■ 

-Ds 

(Mp, 

+ 

MP6) 

(3.14) 


This  transfer  function  is  depicted  in  the  signal  flow  graph  of  Figure  6  on  page  16. 
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Figure  6.      Pitch  Transfer  Function  Signal  Flow  Graph 

The  aerodynamic  damping  term,  D,  is  calculated  for  each  of  the  three  design  points. 
The  denominator  of  the  6(s)l5(s)  transfer  function  is  the  characteristic  equation  of  that 
system.    It  is  in  the  quadratic  form  of  s7  +  s2((on  +  co*.   Equating  coefficients  yields 


-  D  =  2{co„ 


(3.15) 


and 


{Mpa.  +  MP6) 


(3.16) 


Using  the  missile  parameters  obtained  in  the  previous  chapter  and  the  assumed  value 
of  0.5  for  (,  the  damping  term,  D,  is  calculated.  The  values  for  D  are  given  in  Table  5 
on  page  17.   [Ref.  5:  pp.  106-107] 


■ 
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Table  5.     AERODYNAMIC     DAMPING     COEFFI 
CI  EN  IS 


DESIGN  POINT 

DAMPING  COEFFICIENT 

1 

-22.4 

2 

-11.9 

3 

-9.1 

The  other  missile  dynamics  transfer  function  needed  is  the  lateral  acceleration 
transfer  function,  7]M(s)l3{s).  Combining  the  signal  How  graph  of  Figure  6  with 
Equation  (3.5)  results  in  the  signal  How  graph  of  Figure  7. 


Figure  7.      Acceleration  Transfer  Function  Signal  Flow  Graph 


The  transfer  function  associated  with  this  signal  flow  graph  is  given  as 


d{s) 


MP6 

Vl*  +  FLt) 

I 

III 

2 
S 

-  sD  ■ 

(MPx  + 

MpJ 

+ 


La. 

m 


(3-17) 
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The  transfer  functions  for  the  missile  dynamics  are  thus  obtained.  They  will  be  used 
in  the  following  section  to  determine  the  autopilot  transfer  functions. 

C.  AUTOPILOT  TRANSFER  FUNCTIONS 

By  inserting  the  two  missile  dynamics  transfer  functions  into  each  of  the  two 
autopilot  configuration  block  diagrams  presented  in  Chapter  One,  and  performing  a  few 
basic  block  diagram  manipulations,  the  closed  loop  autopilot  transfer  functions  are  ob- 
tained. In  both  cases,  the  desired  autopilot  transfer  function  is  rjM(s)lfjc(s).  For  Autopilot 
Configuration  One  the  transfer  function  is 


Yl  v(s)  S(s) 

J^TL  = : — : •  (3.18) 

Ids)         2  ,      6{s)  ,  d{s) 

s  +  s(KF+  KBR  -t—  )  +  KSTKBR  -r—  +  KA 

Expanding  the  missile  dynamics  transfer  functions  in  the  autopilot  transfer  function  and 
performing  extensive  algebraic  manipulations  yields 

G,(5)  =  — tt  =  ~ } 2 O-19) 

;/d5J  5     +<2f,i5     +  d,,S    +  d-),S  +  d-x 


'21-*  ^  "31 


where 


nM  = 


FL6 


01  -     m 


nn=-D 


m 


MP6   (FL,  +  FL6)       FL6   (MPx  +  MPd) 

n-y-i  = — : 


m  in 


d0l  =  KF-  D, 


MP5  FL6  (A/„,  +  MPS) 


'\\  —  1XBR       j       i"  ^A     m     ~  ^¥u  ~  j 


.     _  K    K     Mn      K   (^  +  Mm)      „  n  Fu 

"21  -  ^ST^BR        7  KF  7  _  liAU~nT~ 


and 


IS 


,        „  ,  MH   {Fu  +  FL<5)       FL6   {MPy  +  MP6) 

The  transfer  function  for  Autopilot  Configuration  Two  is  obtained  in  a  similar 
manner.    For  the  second  configuration 


(KyS+K,) 


tlu(s)  2  '      6(s 


where 


vh*. 

w02  —  A2     /7/       ' 


.      FL6  „    r,    F> 


»n  =  A.',      ,,,     —  A'-,/) 


Z.<5 


(3.20) 


*CW  2  0($)  >7V/(S) 

This  expands  into 

r<\        y]U[s)  n02s3  +  "\2s2  +  "22s  +  "32  ,,  0f, 

G2(s)  = — TT  =  _4 5 5 ^3-21^ 

few        s   +  d02s ■  +dns +  dns  +  d32 


2      /vi    m  2       m     ' 

„,  MP6    (Flz  +  Fls)       fl6    (MPx  +  MP6)  %  i^ 

«22-A2l       y  m  -     m  7  )-^l^     m      , 

«32~AA      j  ,n         -    m  j  ), 

FL6 

^02  =  kf  -  D  +  KAK2  —j^  , 

MP6       Mp*  +  Mpfi        ,  ,      FL6  r   FLfi 

"12  -  kbr  ~  ~i  KFD  -  KAK2D  —ppj-  +  KAKl  -fij-  , 

,        v  vt  MP6   (FL,  +  FL6)       FL6   (MPy  +  Mps)  Fu  (A/Pg  +  Afw) 

"22-A^^2l       J  m  _     m  J  )-  ^AK\U~JT:  KF  J 

and 
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A  »    *>  <   M™     {FL*  +  FL8)  FL6     {MP*  +  MP6)    ^ 

"32 —  *\4ah      j  m  m  j  >  ' 

These  closed  loop  autopilot  transfer  functions  will  be  used  in  the  next  section  for 
system  pole  placement  and  solution  of  the  feedback  and  feed  forward  gains. 

D.     POLE  PLACEMENT  AND  GAIN  DETERMINATION 

As  stated  in  Chapter  One  of  this  research  paper,  the  desired  autopilot  response  for 
both  configurations  is  con  =  10  with  a  damping  coefficient  £  =  0.5  .  These  design  spec- 
ifications dictate  a  complex  pair  of  poles  in  the  characteristic  equation  whose  real  part 
is  -5.0.  Since  the  autopilots'  characteristic  equations  are  fourth  order,  two  more  poles 
are  required  to  put  them  in  the  form 

CE  =  (s2  +  s2£con  +  (o2n)(s  +  Pl){s  +  p2) .  (3.22) 

In  order  for  the  complex  pair  placed  by  design  specifications  to  be  dominant,  the  other 
two  poles  must  lie  to  the  left  of  them  in  the  s-plane.  The  two  additional  poles  are  se- 
lected for  simplicity  as  real  and  equal  and  placed  at  -20.0.    The  characteristic  equation 

for  both  autopilot  configurations  is 

CE={s+  20)(s  +  20)02  +  10s  +  100) 

=  /  +  50s3  +  900s2  +  SOOOs  +  40000  .  (3.23) 

The  denominators  of  the  closed  loop  transfer  functions,  which  are  given  by 
Equations  (3.19)  and  (3.21),  are  the  autopilot  characteristic  equations.  The  transfer 
function  gains  are  solved  for  by  first  equating  the  coefficients  of  characteristic  equation 
of  above  to  those  of  Equations  (3.19)  and  (3.21),  then  solving  simultaneous  equations 
for  the  gains.  The  calculated  gains  for  Autopilot  Configuration  One  and  Autopilot 
Configuration  Two  are  tabulated  in  Table  6  and  Table  7,  respectively.  The  determi- 
nation of  the  autopilot  gains  makes  it  possible  to  simulate  the  response  of  the  missile. 
These  simulations  will  be  conducted  in  the  next  chapter.   [Refs.  6,7] 
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Table  6.     AUTOPILOT  GAINS  FOR  CONFIGURATION  ONE 


GAIN 

DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 
POINT  3 

Ka 

-0.094 

-0.S45 

-3.06 

Kf 

27.6 

38.1 

40.9 

KST 

26.96 

8.9S 

10.18 

KBR 

0.752 

-4.44 

-12.37 

Table  7.     AUTOPILOT  GAINS  FOR  CONFIGURATION  TWO 


GAIN 

DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 
POINT  3 

KAKX 

-0.094 

-0.845 

-3.06 

KAK2 

0.004 

-0.057 

-0.352 

Kf 

20.1 

41.05 

51.17 

Kbr 

-0.513 

-7.62 

-12.83 

It  is  interesting  to  note  that  the  values  of  A',  and  A',  for  Configuration  Two,  depend 
on  the  selected  value  of  A'^.  For  purposes  of  simulation  in  the  next  chapter,  A^  is  chosen 
as  unity.  This  results  in  KAKX  being  equal  to  A',  and  KAK2  being  equal  to  K2  for  the  sim- 
ulation. 
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IV.     AUTOPILOT  SIMULATION  AND  ANALYSIS 

A.     AUTOPILOT  SIMULATION 

In  order  to  compare  the  responses  of  the  two  autopilot  configurations,  a  simulation 
program  was  developed  using  IBM's  Dynamic  Simulation  Language.   [Ref.  8] 

.  The  program  simulates  the  step  response  of  the  autopilot  transfer  functions  for  both 
configurations,  at  all  three  design  points.  The  simulation  was  developed  from  Equations 
(3.19)  and  (3.21),  which  describe  the  two  autopilot  transfer  functions,  and  uses  the 
feedback  and  feed  forward  gains  tabulated  in  Table  6  and  Table  7  on  page  21.  The 
program  also  incorporates  a  saturating  limiter  which  limits  the  deflection  of  the  control 
fins  to  ±  20°.  A  copy  of  the  program  listing  is  at  Appendix  E. 

The  program  was  run  on  the  IBM  3360  mainframe  computer  system  at  the  Naval 
Postgraduate  School  utilizing  a  Tektronix  618  Monitor  and  Tektronix  4631  Hard  Copy 
Unit  for  graphic  output. 

This  chapter  contains  the  graphical  results  of  these  simulations.  Figure  8  on  page 
23  through  Figure  10  on  page  25  depict  the  step  response  of  Autopilot  Configuration 
One  at  each  of  the  three  design  points.  Figure  1 1  on  page  26  shows  the  parametric 
display  of  the  step  response  of  Autopilot  Configuration  One  at  all  three  design  points. 
Figure  12  on  page  27  through  Figure  14  on  page  29  depict  the  step  response  of 
Autopilot  Configuration  Two  at  each  of  the  three  design  points.  Figure  15  on  page  30 
displays  parametrically  the  step  response  of  Autopilot  Configuration  Two  at  all  three 
design  points.  The  parametric  comparisons  of  the  two  autopilot  configurations  at  each 
design  point  are  shown  in  Figure  16  on  page  31  through  Figure  IS  on  page  33.  The 
next  section  contains  discussion  and  analysis  of  the  graphical  output. 


*)-> 


Figure  8.      Step  Response  of  Autopilot  Configuration  One  at  Design  Point  One 
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Figure  9.      Step  Response  of  Autopilot  Configuration  One  at  Design  Point  Two 
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Figure  10.      Step  Response  of  Autopilot  Configuration  One  at  Design  Point  Three 
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Figure  11.      Parametric  Comparison  of  the  Step  Response  of  Autopilot  Configuration 
One  at  Three  Design  Points 
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Figure  12.      Step  Response  of  Autopilot  Configuration  Two  at  Design  Point  One 
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Figure  13.      Step  Response  of  Autopilot  Configuration  T>\o  at  Design  Point  Two 
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Figure  14.      Step  Response  of  Autopilot  Configuration  Two  at  Design  Point  Three 
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Figure  15.      Parametric  Comparison  of  the  Step  Response  of  Autopilot  Configuration 
Two  at  Three  Design  Points 
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Figure  16.      Parametric  Comparison  of  the  Step  Response  of  Autopilot  Configura- 
tions One  and  Two  at  Design  Point  One 
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Figure  17.       Parametric  Comparison  of  the  Step  Response  of  Autopilot  Configura- 
tions One  and  T»o  at  Design  Point  Two 
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Figure  18.      Parametric  Comparison  of  the  Step  Response  of  Autopilot  Configura- 
tions One  and  Two  at  Design  Point  Three 
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B.     ANALYSIS  OF  SIMULATION  RESULTS 

In  the  analysis  of  the  simulation  results,  there  are  two  characteristics  of  the  step  re- 
sponse which  are  of  interest.    They  are  the  rise  time,  t,  ,  and  the  settling  time,  t.  . 

The  rise  time  is  simply  defined  as  the  time  required  for  the  systems  output  response 
to  go  from  10  percent  to  90  percent  of  its  final  value.  This  can  be  measured  directly  off 
of  the  graphical  simulation  results.   [Ref  9:  p.  40] 

•  The  settling  time,  ts  ,  is  defined  as  four  time  constants,  or 

<s  =  T7--  (4-D 

The  value  of  the  maximum  overshoot,  Mp[ ,  and  the  time  at  which  the  maximum  over- 
shoot occurs.  ip  ,  can  be  measured  directly  off  the  graphical  output.  The  maximum 
overshoot  is  defined  as 


A/„=l+exp(      -JH-  )  (4.2) 

\  l  ~  s 


and  the  time  at  which  it  occurs 


t  = £==-.  (4.3) 


&>«\  1 


Measurements  of  Mp,  and  tp  are  made  from  the  graphical  output  of  Figure  8  on  page 
23  through  Figure  10  on  page  25  for  Configuration  One  and  from  Figure  12  on  page 
27through  Figure  14  on  page  29  for  Configuration  Two.  Once  Mpt  and  ip  are  measured, 
£  and  co„  are  obtained  by  solving  simultaneously  Equations  (4.2)  and  (4.3).  The  values 
for  I  and  co„  are  then  used  in  Equation  (4.1)  to  obtain  the  settling  time,  tr  The  response 
characteristics  Mpt,  tp,  C,  con  ,  /,  and  tr  are  tabulated  in  Table  8  and  Table  9  on  page  35. 
[Ref  5:  pp.  10S-109] 
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Table  8.     RESPONSE    CHARACTERISTICS    FOR   AUTOPILOT 
CONFIGURATION  ONE 


DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 

POINT  3 

Mpl 

1.12 

1.13 

1.07 

h 

0.51  sec 

0.51  sec 

0.57  sec 

0.56 

0.54 

0.65 

Q)n 

7.47  rad'sec 

7.36  radsec 

7.24  rad'sec 

ls 

0.96  sec 

1.01  sec 

0.85  sec 

h 

0.21  sec 

0.19  sec 

0.23  sec 

Table  9.     RESPONSE   CHARACTERISTICS    FOR    AUTOPILOT 
CONFIGURATION  TWO 


DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 
POINT  3 

Mpt 

1.13 

1.19 

1.24 

h 

0.63  sec 

0.52  sec 

0.43  sec 

C 

0.54 

0.47 

0.41 

con 

5.95  rad  sec 

6.84  rad  sec 

8.03  rad  sec 

<s 

1.24  sec 

1.24  sec 

1.20  sec 

*r 

0.23  sec 

0.19  sec 

0.14  sec 
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V.     CONCLUSIONS 

The  step  responses  of  the  two  autopilot  configurations  were  compared  at  three  dif- 
ferent design  points.  Both  configurations  met  the  desired  response  specifications  of 
C^0.5  and  cor^\0.  The  actual  values  of  £  and  co„  are  given  in  Table  8  and  Table  9  on 
page  35  for  Configurations  One  and  Two  respectively. 

At  Design  Point  One,  Configuration  One's  observed  rise  time  is  nine  percent  less 
than  that  of  Configuration  Two,  and  Configuration  One's  observed  settling  time  is 
twenty-three  percent  less  than  that  of  Configuration  Two.  The  two  configurations  are 
compared  graphically  in  Figure  16  on  page  31. 

At  Design  Point  Two.  Configuration  One's  observed  rise  time  is  equal  to  that  of 
Configuration  Two.  and  Configuration  One's  observed  settling  time  is  nineteen  percent 
less  than  that  of  Configuration  Two.  The  graphical  comparison  is  shown  in  Figure  17 
on  page  32. 

At  Design  Point  Three,  Configuration  One's  observed  rise  time  is  sixty-four  percent 
more  than  that  of  Configuration  Two,  and  Configuration  One's  observed  settling  time 
is  thirty  percent  less  than  that  of  Configuration  Two.  This  comparison  is  shown  graph- 
ically in  Figure  18  on  page  33. 

The  apparent  trend  is  that  Configuration  One's  rise  time  increases  slightly  as  the 
flight  of  the  missile  continues  but  that  its  settling  time  decreases  slightly.  This  trend  is 
graphically  depicted  in  Figure  11  on  page  26.  Configuration  Two's  rise  time  decreases 
significantly  as  the  flight  of  the  missile  continues  but  its  settling  time  remains  nearly 
constant.    This  trend  is  shown  graphically  in  Figure  15  on  page  30. 

The  shorter  rise  time,  later  in  flight,  of  Configuration  Two.  is  viewed  as  a  significant 
performance  advantage  over  Configuration  One.  It  is  clear  that  a  shorter  rise  time 
means  a  faster  response.  This  faster  response  is  desirable  late  in  flight  because  it  is 
during  this  later  part  of  the  missile's  flight  that  the  missile  would  be  in  its  final  phase  of 
closing  on  and  tracking  a  maneuvering  target.  For  this  reason,  faster  response,  later  in 
flight,  is  a  significant  advantage. 
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APPENDIX  A.     ORIGINAL  WIND  TUNNEL  DATA 


Table   10.     AXIAL  FORCE  COEFFICIENTS  (POWER  ON) 

ALT  (KM) 

/ 1           \ 

ALPHA 

MACH  2.60 

MACH  3.20 

MACH  3.80 

(km) 

(degrees) 

0 

0.324 

0.2SS 

0.25S 

3.0 

0.326 

0.289 

0.260 

0 

5.0 

0.328 

0.292 

0.263 

10.0 

0.338 

0.305 

0.278 

20.0 

0.358 

0.336 

0.322 

30.0 

0.248 

0.323 

0.422 

o 

0.334 

0.295 

0.263 

3.0 

0.335 

0.296 

0.265 

8 

5.0 

0.33S 

0.299 

0.26S 

10.0 

0.348 

0.312 

0.283 

20.0 

0.367 

0.343 

0.328 

30.0 

0.257 

0.330 

0.427 

0 

0.353 

0.311 

0.278 

3.0 

0.355 

0.313 

0.279 

16 

5.0 

0.357 

0.316 

0.282 

10.0 

0.367 

0.329 

0.297 

20.0 

0.387 

0.360 

0.342 

30.0 

0.276 

0.347 

0.441 

0 

0.381 

0.336 

0.299 

3.0 

0.383 

0.337 

0.301 

24 

5.0 

0.385 

0.340 

0.304 

10.0 

0.395 

0.353 

0.319 

20.0 

0.415 

0.3S4 

0.363 

30.0 

0.304 

0.371 

0.462 
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Table   11.     AXIAL  FORCE  COEFFICIENTS  (POWER  OFF) 

ALT 

ALPHA 

MACH  2.60 

MACH  3.20 

MACH  3.80 

(km) 

(degrees) 

0 

0.392 

0.342 

0.303 

3.0 

0.394 

0.344 

o.3o5 

0 

5.0 

0.396 

0.347 

0.308 

10.0 

0.407 

0.360 

0.323 

20.0 

0.429 

0.391 

0.367 

30.0 

0.321 

0.378 

0.466 

0 

0.401 

0.349 

0.30S 

3.0 

0.403 

0.351 

0.310 

8 

5.0 

0.406 

0.354 

0.313 

10.0 

0.416 

0.367 

0.328 

20.0 

0.43S 

0.39S 

0.372 

30.0 

0.331 

0.385 

0.471 

0 

0.421 

0.366 

0.322 

3.0 

0.422 

0.368 

0.324 

16 

5.0 

0.425 

0.370 

0.327 

10.0 

0.436 

0.383 

0.342 

20.0 

0.458 

0.414 

0.387 

30.0 

0.350 

0.401 

0.486 

0 

0.449 

0.390 

0.344 

3.0 

0.451 

0.392 

0.345 

24 

5.0 

0.453 

0.395 

0.348 

10.0 

0.464 

0.408 

0.363 

20.0 

0.4S6 

0.438 

0.40S 

30.0 

0.378 

0.425 

0.507 
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Table  12.     NORMAL  FORCE  COEFFICIENTS 


DELTA 

(degrees) 

ALPHA 

(degrees) 

MACH  2.60 

MACH  3.20 

MACH  3.80 

0 

-1.104 

-0.999 

-0.525 

3.0 

-0.459 

-0.392 

-0.015 

-30 

5.0 

-0.047 

0.013 

0.353 

10.0 

1.336 

1.363 

1.594 

20.0 

5.004 

4.877 

4.739 

30.0 

9.420 

8.992 

8.391 

0 

-0.802 

-0.660 

-0.423 

3.0 

-0.138 

-0.041 

O.102 

-20 

5.0 

0.291 

0.368 

0.479 

10.0 

1.699 

1.707 

1.743 

20.0 

5.358 

5.163 

4.926 

30.0 

9.758 

9.259 

8.604 

0 

-0.3S5 

-0.305 

-0.236 

3.0 

0.272 

0.294 

0.300 

-10 

5.0 

0.689 

0.692 

0.683 

10.0 

2.076 

2.008 

1.959 

20.0 

5.723 

5.451 

5.159 

30.0 

10.156 

9.577 

8.839 

0 

0.000 

o.ooo 

0.000 

3.0 

0.652 

0.594 

0.539 

0 

5.0 

1.067 

0.989 

0.923 

10.0 

2.457 

2.309 

2.201 

20.0 

6.142 

5.7SS 

5.391 

30.0 

10.025 

10.004 

9.043 

0 

0.385 

0.305 

0.236 

3.0 

1.046 

0.908 

0.770 

10 

5.0 

1.469 

1.312 

1.151 

10.0 

2.888 

2.660 

2.471 

20.0 

6.585 

6.225 

5.572 

30.0 

11.014 

10.421 

9.175 

0 

0.802 

0.660 

0.423 

3.0 

1.460 

1.286 

0.946 

20 

5.0 

1.S79 

1.709 

1.319 

10.0 

3.271 

3.063 

2.563 

20.0 

6.911 

6.574 

5.667 

30.0 

11.264 

10.731 

9.208 

0 

1.104 

0.999 

0.525 

3.0 

1.747 

1.606 

1.033 

30 

5.0 

2.154 

2.009 

1.396 

10.0 

3.513 

3.342 

2.613 

20.0 

7.072 

6.804 

5.659 

30.0 

11.430 

10.889 

9.141 
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Table   13.     PITCHING  MOMENT  COEFFICIENTS 

DELTA 

ALPHA 

MACH  2.60 

MACH  3.20 

MACH  3.80 

(desrees) 

(degrees) 

0 

7.313 

6.609 

3.437 

3.0 

7.308 

6.480 

3.601 

-30 

5.0 

7.016 

6.14S 

3.451 

10.0 

6.632 

5.432 

3.254 

20.0 

6.353 

4.212 

2.740 

30.0 

2.980 

0.503 

1.222 

0 

5.307 

4.366 

2.797 

3.0 

5.185 

4.161 

2.828 

-20 

5.0 

4.782 

3.802 

2.616 

10.0 

4.227 

3.152 

2.272 

20.0 

4.013 

2.317 

1.508 

30.0 

0.745 

-1.259 

-0.188 

0 

2.550 

2.017 

1.560 

3.0 

2.469 

1.944 

1.524 

-10 

5.0 

2.146 

1.662 

1.271 

10.0 

1.735 

1.160 

0.839 

20. 0 

1.592 

0.414 

-0.034 

30.0 

-1.S89 

-3.362 

-1.736 

0 

o.ooo 

0.000 

0.000 

3.0 

-0.043 

-0.034 

-0.054 

0 

5.0 

-0.355 

-0.303 

-0.315 

10.0 

-0.780 

-0.818 

-0.749 

20.0 

-1.167 

-1.805 

-1.553 

30.0 

-5.033 

-6.229 

-3.093 

0 

-2.550 

-2.017 

-1.560 

3.0 

-2.655 

-2.115 

-1.586 

10 

5.0 

-3.023 

-2.441 

-1.824 

10.0 

-3.646 

-3.150 

-2.184 

20.0 

-4.113 

-4.706 

-2.764 

30.0 

-7.570 

-8.942 

-3.960 

0 

-5.307 

-4.366 

-2.797 

3.0 

-5.398 

-4.618 

-2.750 

20 

5.0 

-5.735 

-5.067 

-2.935 

10.0 

-6.182 

-5.816 

-3.149 

2o.O 

-6.269 

-7.017 

-3.388 

30.0 

-9.224 

-10.995 

-4.178 

0 

-7.313 

-6.609 

-3.  4  /  j 

3.0 

-7.297 

-6.732 

-3.326 

30 

5.0 

-7.556 

-7.054 

-3.442 

10.0 

-7.784 

-7.664 

-3.479 

20.0 

-7.338 

-S.535 

-3.339 

30.0 

-10.323 

-12.038 

-3.732 
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APPENDIX  B.     INTERPOLATED  WIND  TUNNEL  DATA 


Table  14.     AXIAL  FORCE  COEFFICIENTS  (POWER  ON) 


ALPHA 

DESIGN 

DESIGN 

DESIGN 

(degrees) 

POINT  1 

POINT  2 

POINT  3 

0 

0.278 

0.280 

0.319 

3.0 

0.2  SO 

0.2S1 

0.320 

5.0 

0.283 

0.284 

0.323 

10.0 

0.296 

0.299 

0.336 

20.0 

0.335 

0.341 

0.367 

30.0 

(J.373 

0.4  IS 

0.349 

Table  15.     AXIAL  FORCE  COEFFICIENTS  (POWER  OFF) 


ALPHA 

DESIGN 

DESIGN 

DESIGN 

(degrees) 

POINT  I 

POINT  2 

POINT  3 

0 

0.310 

0.326 

0.374 

3.0 

0.329 

0.328 

0.376 

5.0 

0.332 

0.331 

0.379 

10.0 

0.347 

0.346 

0.392 

20.0 

0.383 

0.387 

0.421 

30.0 

0.421 

0.464 

0.404 
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Table  16.     NORMAL  FORCE  COEFFICIENTS 


DELTA 

(degrees) 

ALPHA 

(decrees) 

DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 
POINT  3 

-30 

0 

3.0 

5.0 

10.0 

20.0 

30. 0 

-0.778 
-0.216 
0.332 
1.417 
4.813 
8.712 

-0.620 
-0.090 
0.2S5 
1.548 
4.767 
8.511 

-1.004 
-0.395 
0.010 
1.362 
4.883 
9.013 

-20 

0 

3.0 
5.0 
10.0 

20.0 
30.0 

-0.549 
0.026 
0.420 
1.724 
5.052 
8.953 

-0.470 
0.073 
0.456 
1.736 

4.973 

8.735 

-0.667 

-0.046 

0.364 

1.707 

5.173 
9.284 

-10 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

-0.273 
0.297 
0.6SS 
1.985 
5.315 
9.233 

-0.250 
0.299 
0.6S5 
1.969 
5.217 
8.987 

-0.309 

0.293 

0.692 

2.011 

5.465 

9.606 

0 

0 

3.0 
5.0 
10.0 
20.0 

30.0 

0.000 
0.568 
0.958 
2.259 
5.603 
9.556 

0.000 
0.550 
0.936 
2.223 

5.470 
9.235 

0.000 

0.597 
0.993 
2.316 
5.S06 
10.035 

10 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

0.273 
0.S44 

1.237 
2.547 
5.920 
9.839 

0.250 
0.798 
1.183 
2.466 
5.703 
9.424 

0.309 
0.915 
1.320 
2.671 
6.243 
10.451 

20 

0 

3.0 
5.0 
10.0 
20.0 
30.0 

0.549 
1.127 
1.527 
2.830 
6.151 
10.020 

0.470 
1.014 
1.397 
2.663 
5.848 
9.513 

0.66S 
1.295 
1.718 

3.073 
6.591 
10.758 

30 

0 

3.0 
5.0 
10.0 
20.0 
30.0 

0.77S 
1.339 
1.723 
3.002 
6.26S 
10.073 

0.620 
1.148 
1.519 
2.759 
5.88S 
9.491 

1.004 
1.613 
2.016 
3.351 
6.817 
10.916 
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Table  17.     PITCHING  MOMENT  COEFFICIENTS 

DELTA 

ALPHA 

DESIGN 

DESIGN- 

DESIGN 

(degrees) 

(degrees) 

POINT  1 

POINT  2 

POINT  3 

0 

5.145 

4.1()() 

6.644 

3.0 

5.136 

4.177 

6.521 

-30 

5.0 

4.889 

3.990 

6.191 

10.0 

4.416 

3.690 

5.492 

20.0 

3.525 

3.036 

4.319 

30.0 

0.839 

1.078 

0.627 

0 

3.634 

3.111 

4.413 

3.0 

3.539 

3.095 

4.212 

-20 

5.0 

3.24S 

2.853 

3.851 

10.0 

2.741 

2.448 

3.206 

20.0 

1.939 

1.670 

2.402 

30.0 

-0.759 

-0.402 

-1.159 

0 

1.804 

1.651 

2.044 

3.0 

1.74S 

1.608 

1.970 

-10 

5.0 
10.0 

1.4S0 
1.010 

1.349 
0.903 

1.686 
1.189 

20.0 

0.205 

0.056 

0.473 

30.0 

-2.603 

-2.061 

-3.2SS 

0 

0.000 

0.000 

0.000 

3.0 

-0.043 

-0.050 

-0.034 

0 

5.0 

-0.309 

-0.313 

-0.306 

10.0 

-0.786 

-0.763 

-0.816 

20.0 

-1.6S7 

-1.603 

-1.773 

30.0 

-4.765 

-3.720 

-6.169 

0 

-1.804 

-1.651 

-2.044 

3.0 

-1.868 

-1.692 

-2.142 

10 

5.0 

-2.153 

-1.947 

-2.470 

10.0 

-2.669 

-2.377 

-3.175 

20.0 

-3.800 

-3.152 

-4.676 

30.0 

-6.617 

-4.956 

-8.873 

0 

-3.634 

-3.111 

-4.413 

3.0 

-3.746 

-3.124 

-4.657 

20 

5.0 

-4.072 

-3.361 

-5.100 

10.0 

-4.571 

-3.680 

-5.834 

20.0 

-5.323 

-4.050 

-6.980 

30.0 

-7.814 

-5.541 

-10.906 

0 

-5.145 

-4.100 

-6.644 

3.0 

-5.142 

-4.007 

-6.760 

30 

5.0 

-5.36S 

-4.164 

-7.079 

10.0 

-5.711 

-4.316 

-7.670 

20.0 

-6.110 

-4.378 

-S.475 

30.0 

-8.162 

-5.393 

-11.952 
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APPENDIX  C.     LIFT  AND  DRAG  COEFFICIENT  DATA 


Table   18.     LIFT  COEFFICIENTS  (POWER  ON) 


DELTA 

ALPHA 

DESIGN 

DESIGN 

DESIGN 

(degrees) 

(degrees) 

POINT  1 

POINT  2 

POINT  3 

0 

-0.778 

-0.620 

- 1 .004 

3.0 

-0.230 

-0.105 

-0.411 

-30 

5.0 

0.306 

0.259 

-0.007 

10.0 

1.397 

1.473 

1.283 

20.0 

4.408 

4.363 

4.463 

30.0 

7.358 

7.162 

7.631 

0 

-0.549 

-0.470 

-0.668 

3.0 

0.011 

0.058 

-0.063 

-20 

5.0 

0.394 

0.429 

0.334 

10.0 

1.646 

1.658 

1.623 

20.0 

4.633 

4.556 

4.736 

30.0 

7.567 

7.356 

7.866 

0 

-0.273 

-0.250 

-0.309 

3.0 

0.282 

0.2S3 

0.276 

-10 

5.0 

0.661 

0.658 

0.661 

10.0 

1.903 

1.887 

1.922 

20.0 

4.S80 

4.786 

5.010 

30.0 

7.810 

7.574 

8.145 

0 

0.000 

0.000 

0.000 

3.0 

0.553 

0.535 

0.579 

0 

5.0 

0.930 

0.908 

0.961 

10.0 

2.173 

2.137 

2  ">")~) 

20.0 

5.151 

5.023 

5. 3  30 

30.0 

S.0S9 

7.789 

S.516 

0 

0.273 

0.250 

0.309 

3.0 

O.S28 

0.7S2 

0.897 

10 

5.0 

1.208 

1.154 

1.287 

10.0 

2.457 

2.377 

~>  S7^ 

20.0 

5.448 

5.242 

5.741 

30.0 

8.334 

7.952 

S.876 

0 

0.549 

0.470 

0.66S 

3.0 

1.111 

0.998 

1.276 

20 

5.0 

1.497 

1.367 

1.683 

10.0 

2.736 

2.571 

2.968 

20.0 

5.665 

5.379 

6.068 

30.0 

8.491 

8.030 

9.142 

0 

0.778 

0.620 

1.004 

3.0 

1.323 

1.132 

1.594 

30 

5.0 

1.692 

1.488 

1.980 

10.0 

2.905 

2.665 

3.242 

20.0 

5.775 

5.416 

6.280 

30.0 

8.537 

8.010 

9.279 
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Table  19.     LIFT  COEFFICIENTS  (POWER  OFF) 


DELTA 

ALPHA 

DESIGN 

DESIGN 

DESIGN- 

(degrees) 

(degrees) 

POINT  1 

POINT  2 

POINT  3 

0 

-0.77S 

-0.620 

-1.004 

3.0 

-n.233 

-0.107 

-0.414 

-30 

5.0 

0.303 

0.255 

-0.023 

10.0 

1.3SS 

1.464 

1.273 

20.0 

4.392 

4.347 

4.445 

30.0 

7.334 

7.139 

7.603 

0 

-0.549 

-0.470 

-0.668 

3.0 

0.009 

0.056 

-0.066 

-20 

5.0 

0.3S9 

0.425 

0.330 

10.0 

1.638 

1.650 

1.613 

20.0 

4.616 

4.541 

4.717 

30.0 

7.543 

7.333 

7.838 

0 

-0.273 

-0.250 

-0.309 

3.0 

0.279 

0.281 

0.273 

-10 

5.0 

0.056 

0.654 

0.656 

10.0 

1.895 

1.879 

1.912 

20.0 

4.S63 

4.770 

4.991 

30.0 

7.786 

7.551 

8.117 

0 

o.ooo 

0.000 

0.000 

3.0 

0.550 

0.532 

0.577 

0 

5.0 

0.925 

0.904 

0.956 

10.0 

2.164 

2.129 

2.213 

20.0 

5.134 

5.00S 

5.312 

30.0 

8.065 

7.766 

8.4S9 

0 

0.273 

0.250 

0.309 

3.0 

0.8265 

0.780 

0.894 

10 

5.0 

1.203 

1.150 

1.282 

10.0 

2.44S 

2.36S 

2.562 

20.0 

5.432 

5.227 

5.723 

30.0 

8.310 

7.929 

8.849 

0 

0.549 

0.4  7o 

0.668 

3.0 

1.108 

0.995 

1.274 

20 

5.0 

1.492 

1.363 

1.67S 

10.0 

2.727 

2.562 

2.958 

20.0 

5.649 

5.363 

6.050 

30.0 

S.467 

8.0075 

9.115 

0 

0.778 

0.620 

1 .004 

3.0 

1.327 

1.129 

1.591 

30 

5.0 

1.688 

1.484 

1.975 

10.0 

2.896 

2.657 

3.232 

20.0 

5.759 

5.401 

6.262 

30.0 

8.513 

7.987 

9.252 
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Table  20.     DRAG  COEFF1 

[CIENTS  (POWER  ON) 

DELTA 

ALPHA 

DESIGN 

DESIGN 

DESIGN 

( decrees) 

(degrees) 

POINT  1 

POINT  2 

POINT  3 

0 

0.27S 

0.2  SO 

0.319 

3.0 

0.268 

0.276 

0.299 

-30 

5.0 

0.311 

0.308 

0.323 

10.0 

0.556 

0.563 

0.567 

20.0 

1.196 

1.951 

2.015 

30.0 

4.679 

4.618 

4.809 

0 

0.278 

0.280 

0.319 
0.317 
0.353 
2.114 

4.944 

3.0 

0.281 

0.284 

-20 

5.0 

0.319 

0.323 

10.0 

0.591 

0.596 

20.0 

2.043 

2.021 

30.0 

4.800 

4.730 

0 

0.278 

0.280 

0.319 

3.0 

0.295 

0.296 

0.335 

-10 

5.0 

0.342 

0.340 

0.382 

10.0 

0.636 

0.636 

0.680 

20.0 

2.133 

2.105 

2.214 

30.0 

4.940 

4.856 

5.105 

0 

0.278 

0.280 

0.319 

3.0 

0.309 

0.309 

0.351 

0 

5.0 

0.365 

0.364 

0.408 

10.0 

0.684 

0.680 

0.733 

20.0 

2.231 

2.191 

2.331 

30.0 

5.101 

4.980 

5.320 

0 

0.278 

0.280 

0.319 

3.0 

0.324 

0.322 

0.367 

10 

5.0 

0.390 

0.3S6 

0.437 

10.0 

0.734 

0.723 

0.795 

20.0 

2.340 

2.271 

2.4SO 

30.0 

5.243 

5.074 

5.528 

0 

0.278 

0.280 

0.319 

3.0 

0.339 

0.334 

0.387 

20 

5.0 

0.415 

0.405 

0.472 

10.0 

0.783 

0.757 

0.865 

20.0 

2.419 

2.321 

2.599 

30.0 

5.333 

5.119 

5.681 

0 

0.278 

0.280 

0.319 

3.0 

0.350 

0.341 

0.404 

30 

5.0 

0.432 

0.415 

0.497 

10.0 

0.813 

0.774 

0.913 

20.0 

2.459 

2.334 

2.676 

30.0 

5.360 

5.108 

5.760 
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Table  21.     DRAG  COEFF1 

CIENTS  (POWER  OFF) 

DELTA 

ALPHA 

DESIGN 

DESIGN 

DESIGN 

(degrees) 

(degrees) 

POINT  1 

POINT  2 

POINT  3 

0 

0.310 

0.326 

0.374 

3.0 

0.317 

0.323 

0.355 

-30 

5.0 

0.360 

0.355 

0.37S 

10.0 

0.597 

0.610 

0.623 

20.0 

2.006 

1.994 

2.066 

30.0 

-4.721 

4.657 

4.856 

0 

0.310 

0.326 

0.374 

3.0 

0.330 

0.331 

0.373 

-20 

5.0 

0.367 

0.369 

0.409 

10.0 

0.641 

0.642 

0.682 

20.0 

2.0S8 

2.065 

2.165 

30.0 

4.841 

4.769 

4.992 

0 

0.310 

0.326 

0.374 

3.0 

0.344 

0.343 

0.391 

-10 

5.0 

0.341 

0.389 

0.43S 

10.0 

0.686 

0.683 

0.735 

20.0 

2.17S 

2.148 

2.265 

30.0 

4.981 

4.895 

5.153 

0 

0.310 

0.326 

0.374 

3.0 

0.358 

0.356 

0.407 

0 

5.0 

0.414 

0.411 

0.464 

10.0 

0.734 

0.727 

0.788 

20.0 

2.276 

2.235 

2.3S1 

30.0 

5.143 

5.019 

5.367 

0 

0.310 

0.326 

0.374 

3.0 

0..373 

0.369 

0.42  3 

10 

5.0 

0.439 

0.433 

0.493 

10.0 

0.7S4 

0.769 

0.850 

20.0 

2.385 

2.314 

2.531 

30.0 

5.284 

5.114 

5.575 

0 

0.310 

0.326 

0.374 

3.0 

0.38S 

0.381 

0.443 

20 

5.0 

0.464 

0.451 

0.527 

10.0 

0.833 

0.803 

0.920 

20.0 

2.464 

2.364 

2.650 

30.0 

5.375 

5.158 

5.729 

0 

0.310 

0.326 

0.0374 

3.0 

0.399 

0.388 

0.460 

30 

5.0 

0.481 

0.462 

0.553 

10.0 

0.863 

0.820 

0.968 

20.0 

2.504 

2.377 

2.727 

30.0 

5.401 

5.147 

5.808 
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APPENDIX  D.     WIND  TUNNEL  GRAPHS 
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Figure  19.      Lift  Coefficient  versus  Angle  of  Attack  at  Design  Point  One 
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Figure  20.      Lift  Coefficient  versus  Fin  Deflection  Angle  at  Design  Point  One 
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Figure  21.      Pitching  Moment  Coefficient  versus  Angle  of  Attack  at  Design  Point 
One 
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Figure  22. 


Pitching  Moment  Coefficient  versus  Fin  Deflection  Angle  at  Design 
Point  One 
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Figure  23.      Lift  Coefficient  versus  Angle  of  Attack  at  Design  Point  T>vo 
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Figure  24.      Lift  Coefficient  versus  Fin  Deflection  Angle  at  Design  Point  Two 
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Figure  25.      Pitching  Moment  Coefficient  versus  Angle  of  Attack  at  Design  Point 
Two 
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Figure  26. 


Pitching  Moment  Coefficient  versus  Fin  Deflection  Angle  at  Design 
Point  Two 
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Figure  27.      Lift  Coefficient  versus  Angle  of  Attack  at  Design  Point  Three 
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Figure  28.      Lift  Coefficient  versus  Fin  Deflection  Angle  at  Design  Point  Three 
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Figure  29.       Pitching  Moment  Coefficient  versus  Angle  of  Attack  at  Design  Point 
Three 
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Figure  30. 


Pitching  Moment  Coefficient  versus  Fin  Deflection  Angle  at  Design 
Point  Three 
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APPENDIX  E.     SIMULATION  COMPUTER  PROGRAM 


f**  THESIS  RESEARCH 


DSL  SIMULATIONS  FOR  TWO  *** 

AUTOPILOT     CONFIGURATIONS  *** 

AT  THREE  DESIGN  POINTS  *** 

JLiJ.A  JL»I>JU 

***                               KENNETH  E.     COCKERHAM  *** 

***                       CAPTAIN,    UNITED   STATES  ARMY  *** 

***                          NAVAL  POSTGRADUATE   SCHOOL  *** 

j-  j-  -♦-  j-  j-  j-  jl .'-  »'-  j-  »'-  »t-  -'*  j-  j-  -*-  ^'-  -*-  -*-  ^  y-  -^  y-  »*-  -'-  *'-  -j-  ~>.  -»-  »*-  -»- »»- »»-  »t- »'-  -♦-  *'-  -'-  J{  «♦-  j-  y-  yf  y-  j-  -»-  *y  jj.  j.  j-  -'-  ju  *'- 


FIXED    INPUT 

*RUN    1 

*****************DESIGN   POINT  ONE******** ************ 


****************M  i  s  S I LE   PARAMETERS******************* 

PARAM   I   =   720.0,    M  =  450.0,    DAMP  =   -22.4,    FLA  =   799276.0,    FLD  =  88808.0 

PARAM  MPA  =   -168736.0,    MPD  =   -191368.0 


********FEEDBACK  GAINS  FOR  CONFIGURATION  ONE********* 
PARAM     KA1   =   -0.094,    KF1   =   27.6,    KBR1   =   0.752,    KST  =   26.96 

********FEEDBACK  GAINS  FOR  CONFIGURATION  TWO********* 

PARAM     KA2   =    1.0,    KF2   =   20.1,    KBR2   =    -0.513 

******FEED  FORWARD  GAINS   FOR  CONFIGURATION  TWO******* 
PARAM      Kl   =    -0.094,    K2   =   0.004 


* 
* 

v.- 
DERI  VAT  I VE 

* 

********* ***AUTOP I LOT  C ONF I GURAT ION  ONE************* 

* 

* 


Rl  =  KA1  *  STEP(0.  0) 

XI  =  Rl    -   X14    -   X12 

X2  =   INTGRL(0.0,X1) 

X3  =  X2    -   Xll    -   X13 
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X4 

= 

INTGRLCO.  0,X3) 

X5 

= 

X4 

DELTA 1 

= 

LIMITC-O.  35,0.  35, X5) 

X6 

= 

DELTA  1  *  MPD/I 

ALFDD1 

= 

X6  +  X7  +  X8 

ALFD1 

= 

INTGRLCO.  0,ALFDD1) 

ALF1 

= 

INTGRL(0.  0,ALFD1) 

X7 

= 

ALFD1  *  DAMP 

X8 

= 

ALF1  *  (MPA+MPD)/I 

X9 

= 

DELTA 1  *  FLD/M 

X10 

= 

ALF1  *  (FLA+FLD)/M 

ACCEL1 

= 

X9  +  X10 

Xll 

= 

KBR1  *  ALFD1 

X12 

= 

KST  *  Xll 

X13 

= 

KF1  *  X4 

X14 

= 

KA1  *  ACCEL1 

* 

* 

* 

* 

*4nfr*********AUT0P I LOT  CONF I GURAT I ON  TWO*************** 


-V 

* 

R2 

= 

KA2  *  STEP(0.  0) 

Yl 

= 

R2  -  Y16 

Y2 

= 

Kl  *  Yl 

Y3 

= 

INTGRLCO.  0,Y2) 

Y4 

= 

K2  *  Yl 

Y5 

= 

Y4  +  Y3 

Y6 

= 

Y5  -  Y14  -  Y15 

Y7 

= 

INTGRLCO.  0,Y6) 

Y8 

= 

Y7 

DELTA2 

= 

LIMITC-O. 35,0. 35, Y8) 

Y9 

= 

DELTA2  *  MPD/I 

ALFDD2 

= 

Y9  +  Y10  +  Yll 

ALFD2 

= 

INTGRLCO. 0,ALFDD2) 

ALF2 

= 

INTGRLCO. 0,ALFD2) 

Y10 

= 

ALF2  *  CMPA+MPD)/I 

Yll 

= 

ALFD2  *  DAMP 

Y12 

= 

ALF2  *  (FLA+FLD)/M 

Y13 

= 

DELTA2  *  FLD/M 

ACCEL2 

= 

Y12  +  Y13 

Y14 

= 

ALFD2  *  KBR2 

Y15 

= 

Y7  *  KF2 

Y16 

= 

KA2  *  ACCEL2 

* 

* 

* 

CONTROL  FINTIM  =2.0,  DELT  =  0.0001 


SAVE  0.001,  ACCEL1,ACCEL2 
END 
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* 

•''RUN   2 

****************DE SIGN  POINT  TWO********************* 

*MrMric4rlck*Mc**#KlSSII£    PARAMETERS******************** 

PARAM  I  =  693,  M  =  370,  DAMP  =  -11.9,  FLA  =  238350.0,  FLD  =  19080.0 
PARAM  MPA  =  -43545.0,  MPD  =  -54431.0 

********FEEDBACK  GAINS  FOR  CONFIGURATION  ONE********* 
PARAM  KA1  =  -0.845,  KF1  =  38.1,  KBR1  =  -4.44,  KST  =  8.98 
********FEEDBACK  GAINS  FOR  CONFIGURATION  TWO********* 
PARAM  KA2  =  1.0,  KF2  =  41.05,  KBR2  =  -7.62 
******F£ED  FORWARD  GAINS  FOR  CONFIGURATION  WO******* 
PARAM   Kl  =  -0.845,  K2  =  -0.057 

END 
* 
* 
* 

*RUN  3 

***************!)£  SIGN  P  0 1  NT  THRE  E  ****^wwr************* 

********m(*****MISSII£   PARAMETERS******************** 

PARAM  I  =  687,  M  =  365,  DAMP  =  -9.1,  FLA  =  99695.0,  FLD  =  10657.0 
PARAM  MPA  =  -26986.0,  MPD  =  -29679.0 

********FEEDBACK  GAINS  FOR  CONFIGURATION  ONE********* 
PARAM  KA1  =  -3.06,  KF1  =  40.9,  KBR1  =  -12.37,  KST  =  10.18 
********FEEDBACK  GAINS  FOR  CONFIGURATION  TWO********* 
PARAM  KA2  =1.0,  KF2  =  51.17,  KBR2  =  -12.83 
******FEED  FORWARD  GAINS  FOR  CONFIGURATION  TWO******* 
PARAM  Kl  =  -3.06,  K2  =  -0.352 

END 
-'- 

GRAPH  (Gl  ,DE=TEK618)  TIME,  ACCEL1(RU=1,LO=-0.  2,SC=0.  2) 
LABEL  (Gl)  STEP  RESPONSE,  CONFIGURATION  1,  DESIGN  POINT  1 

GRAPH  (G2  ,DE=TEK618)  TIME,  ACCEL1(RU=2) 

LABEL  (G2)  STEP  RESPONSE,  CONFIGURATION  1,  DESIGN  POINT  2 

GRAPH  (G3  ,DE=TEK618)  TIME,  ACCEL1(RU=3) 

LABEL  (G3)  STEP  RESPONSE,  CONFIGURATION  1,  DESIGN  POINT  3 

GRAPH  (G4  ,DE=TEK618)  TIME,  ACCEL1(RU=1 , 2 , 3) 

LABEL  (G4)  STEP  RESPONSE,  CONFIGURATION  1,  DESIGN  POINTS  1,  2,  &  3 

GRAPH  (G5  ,DE=TEK618)  TIME,  ACCEL2(RU=1) 

LABEL  (G5)  STEP  RESPONSE,  CONFIGURATION  2,  DESIGN  POINT  1 

GRAPH  (G6  ,DE=TEK618)  TIME,  ACCEL2(RU=2) 

LABEL  (G6)  STEP  RESPONSE,  CONFIGURATION  2,  DESIGN  POINT  2 

GRAPH  (G7  ,DE=TEK618)  TIME,  ACCEL2(RU=3) 

LABEL  (G7)  STEP  RESPONSE,  CONFIGURATION  2,  DESIGN  POINT  3 

GRAPH  (G8  ,DE=TEK618)  TIME,  ACCEL2(RU=1 ,2 , 3) 


62 


LABEL  (G8)  STEP  RESPONSE,  CONFIGURATION  2, DESIGN  POINTS  1,  2,  &  3 

'.V 

GRAPH  (G9  ,DE=TEK618)  TIME,  ACCEL1(RU=1 ,LO=-0. 2 ,SC=0. 2) ,ACCEL2(RU=1) 
LABEL  (G9)  STEP  RESPONSE,  CONFIGURATIONS  1  &  2, DESIGN  POINT  1 

GRAPH  (G10  ,DE=TEK618)  TIME,  ACCEL1(RU=2) , ACCEL2(RU=2) 

LABEL  (G10)  STEP  RESPONSE,  CONFIGURATIONS  1  &  2,  DESIGN  POINT  2 
-v 

GRAPH  (Gil  ,DE=TEK618)  TIME,  ACCEL 1 ( RU=3 ) ,ACCEL2(RU=3) 

LABEL  (Gil)  STEP  RESPONSE,  CONFIGURATIONS  1  6c  2,  DESIGN  POINT  3 

END 
STOP 
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